pacman::p_load(reshape, ggplot2, dplyr, magrittr)
a1 <- "genus.xls" %>% read.delim() %>% melt()
a2 <- "group.txt" %>% read.delim()
a4 <- NULL
for (i in seq_len(nrow(a1))) { a4[i] <- a2[which(a2[, 1] == a1[i, 2]), 2] }
a1[, 4] <- a4


